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Abstract. The classical Hertz entropy is the logarithm of the volume of phase space bounded by the con- 
stant energy surface; its quantum counterpart, the quantum Hertz entropy, is = IniV, i.e. the logarithm 
of the quantum operator A^", specifying the number of states with energy below a given energy eigenstate. 
It has been recently proved that, when an isolated quantum mechanical system is driven out of equilibrium 
by an external driving, the change in the expectation of its quantum Hertz entropy cannot be negative, 
and is null for adiabatic driving. This is in full agreement with the Clausius principle. Here we test the 
behavior of the expectation of the quantum Hertz entropy in the case when two identical XY spin chains 
initially at different temperatures are quenched into a single XY chain. We observed no quantum Hertz 
entropy decrease. This finding further supports the statement that the quantum Hertz entropy is a proper 
entropy for isolated quantum systems. We further quantify how far the quenched chain is from thermal 
equilibrium and the temperature of the closest equilibrium. 

PACS. 05.30.Ch Quantum ensemble theory - 05. 70. -a Thermodynamics - 65.40.Gr Entropy and other 
thermodynamical quantities 



1 Introduction 

The recent tremendous development in the field of nonequi- 
librium quantum fluctuations [HE], has unveiled with an 
unprecedented clarity that many phenomena traditionally 
associated exclusively with macroscopic thermodynamic 
behaviour may manifest themselves even at the micro- 
scopic quantum level. Notably, the Second Law of ther- 
modynamics, in the work- free energy formulation [3^, 

{w)>AF, (1) 

holds down to the quantum level [T. In Eq. ([T]), is the 
work done on a quantum system that is initially in equilib- 
rium with a thermal bath, when the system is perturbed 
by an external time dependent protocol that changes its 
Hamiltonian in time. The brackets indicate average over 
many realizations, and AF is the difference between the 
free energy of a hypothetical equilibrium state (not nec- 
essarily reached by the system) corresponding to the final 
Hamiltonian, and the actual initial free energy of the sys- 
tem. Eq. (IT]) follows straightforwardly from the quantum 
version of the Jarzynski identity [5l|6][7l|8] . 

For a cyclical driving, H{t) = i7(0), AF = 0, Eq. ^ 
says, in accordance with Kelvin postulate, that no energy 
can be extracted by the cyclic variation of a parameter 
from a system in contact with a single bath [4]: 

(w) > Kelvin (cyclic, with bath) (2) 



Given the recent theoretical and experimental advances 
concerning the nonequilibrium dynamics of isolated quan- 
tum systems [9], an interesting question is whether a mi- 
croscopic quantum formulation of the second law in ac- 
cordance with Clausius formulation, is possible as well. 
According to Clausius' formulation the change of entropy 
of a thermally insulated driven system, which begins and 
ends in equilibrium, is non- negative: 

Z\6' > Clausius (no bath) (3) 

Answering this question is not a simple task because it 
amounts to singling out a quantum mechanical quantity 
that behaves as prescribed by the Clausius principle and 
goes over to the usual thermodynamic entropy in the clas- 
sical/thermodynamic limit. Of course von Neumann "en- 
tropy", — Trplogp, proves inadequate in this respect be- 
cause it is invariant under the quantum unitary time evo- 
lution. 

One proposal in the direction of answering the above 
question was put forward by Polkovnikov plQl fTT ] , with the 
introduction of the so-called diagonal entropy. Another 
proposal, put forward by Hal Tasaki and by one of us 
plfni fTI]. uses instead the Hertz entropy |Ttl[T6l [T7] see 
Eq. Q, or in quantum mechanics, its quantum counter- 
part, that is the logarithm of the "quantum number oper- 
ator" N [IMMM^ see Eq. ^ below. In [12,13,14 it was 
shown that the expectation of the quantum Hertz entropy 
behaves in accordance to the Clausius principle. Here we 
scrutinize whether it complies also with another property 
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of thermodynamic entropy, namely whether it increases in 
a scenario when two quantum systems initially at differ- 
ent temperatures are allowed to exchange energy via the 
sudden switch-on of an interaction. This is a scenario that 
has recently attracted considerable attention [I9l[20], and 
that, given the recent advances, e.g., in ultra-cold-atom 
physics [21 , is amenable to experimental investigations. 

In our study the two interacting bodies are two isotropic 
XY spin chains of length A^/2, which are initially at differ- 
ent temperatures, and are suddenly quenched into a single 
isotropic XY spin chain of length N . While the quench dy- 
namics in spin-chains has been thoroughly studied, only 
few studies addressed their thermodynamics [22','23J. 

In Sees. [2] and [3] we review the quantum Hertz entropy 
and present our model, respectively. In Sees. |4] and [5] we 
calculate the initial and final expectation of the quantum 
Hertz entropy. Results concerning the entropy change and 
the deviation of the final state from thermal equilibrium 



where 



and Sec. [7| respectively. Conclu- 



are presented in Sec. [6 
sions are drawn in Sec. 



2 Quantum Hertz entropy 

The microcanonical entropy of a classical system is the so- 
called Hertz entropy p!5 l [T6 l[T7| [T8] presented also by Gibbs 
in his classic book [24], namely 



S = In^/n, 



(4) 



where ^ is the phase space volume enclosed by the hyper- 
surface of constant energy in the system phase space and 
h is Plank's constant. According to semiclassical reason- 
ing [25|, in the quantum limit, the quantity ^/h gets the 
discrete values n = 1,2, . . . The associated quantum 
operator is the number quantum operator TV, whose eigen- 
vectors are the energy eigenvectors, with the integers the 
corresponding eigenvalues. For a driven system, the op- 
erator N is time dependent, and its spectral decompo- 
sition reads (for non degenerate Hamiltonians) N(t) = 
^nnn{t)^ with iln(t) the instantaneous eigenprojectors 
on the eigenspace spanned by the instantaneous eigenvalue 
Enif) of the Hamiltonian H{t). Here it is assumed that the 
energy eigenvalues are non-degenerate and ordered in in- 
creasing fashion: Eiit) < E2{t) < . . . . Accordingly, the 
quantum mechanical operator associated to the Hertz en- 
tropy is 



S{t) =\nN{t). 



(5) 



Under the assumption that (i) the density matrix p(0) 
describing the initial state of the system is diagonal in the 
energy eigenbasis p(0) = ^Pn^n(O), and (ii) the popula- 
tion pn decreases with increasing energy (i.e., Pn < Pm^ if 
En{^) > Em{0)); it has been shown that [T2l[ni [T4]: 



S{t) > S{0) , 



(6) 



^ Depending on the problem at hand, the quantization rule 
may prescribes a shift ^ n + a, which is not relevant in 
this context. 



S{t) = TTS{t)p{t) . 



(7) 



For adiabatic transformations the equal sign holds in Eq. 
([6|. The quantity S{t) hence behaves in accordance to 
the Glausius principle and goes over to the usual ther- 
modynamic entropy in case of large classical systems at 
equilibrium. These facts make it a sound quantum me- 
chanical counterpart for the thermodynamic entropy of a 
thermally insulated system. In the following we shall refer 
to S and S as to the quantum entropy, and the quantum 
entropy operator, respectively. Note however that unless 
the system is in equilibrium at time t, the quantum en- 
tropy S{t) should not be considered as the the thermody- 
namic entropy of the system. The latter is an exclusively 
equilibrium property. 



3 The Model 

Our model consists of two identical isotropic XY spin 
chains, which are quenched at time t = to a single XY 
chain of twice the length, see Fig. [l] At times t < 0, the 
system Hamiltonian is: 



Ho — Hl + Hr , 



(8) 



with: 



N/2 N/2-1 

^^ = "2^-1+2 E [-J-j+i+^H+i]' (9) 

h ^ T 

j=N/2-\-l j^N/2+1 

(10) 



Hq = H^_ + H^ 



Quench at t = 




Hi 

Fig. 1. Schematics of our model 
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Here cr^, j = 1 . . . A^, a = x^y^z^ denotes the PauH ma- 
trices of the j-th spin. At t < the left (right) chain is 
in the Gibbs state of temperature Ti^(^r)^ hence the total 
system density matrix is given by their product: 



-/3lHl 



(11) 



where Z{f3) = Tre~^^^ = Tre~^^^ is the partition func- 
tion, and Pl{r) = {]^bTl(^b))~^ ^ with ks Boltzmann con- 
stant. At time t = an interaction between spin N/2 and 
spin N/2 + 1 is turned on, such that the Hamiltonian is, 
for t > 0: 



N 



J 



N-1 



The Hamiltonians Hl^Hr^Hi all represent isotropic 
XY chains of different lengths. Following the standard 
procedures they can be put in diagonal form by means of 
Jordan- Wigner transformation followed by a sine trans- 
form [26 ^ .27 j . Specifically 



Similarly, for the L-chain: 



N/2 



Nh 



-h — 2 J cos 



kir 



7V/2 + 1 



N/2- 



— N/2 



sm 



kiiT 



N/2+1 



(22) 
(23) 
(24) 



Note that the same Jordan- Wigner operators are used 
for the L-chain and the total chain. In defining all N 
operators are used, while only the first N/2 are employed 
to define the primed operators 6^. We shall denote the 



eigenvectors and eigenvalues of Hl as |/) = 
and Ei^ respectively: 

N/2 

Hl\1)=Ei\1), Ei = Y,^kh 



Likewise for the i?-chain, 



Nh 



(25) 



AT 



where 



Nh 



-h — 2 J cos 



N 



N- 



N- 



1^ 



sm 



kin 
7V + 1 



n 

k=i 



(13) 



(14) 
(15) 

(16) 
(17) 



We shall denote the eigenvectors and eigenvalues of Hi 
as \n) = |ni, 712, . . . ^at) and En^ respectively, where Ui = 
0,1: 



N 



Hi\n) Er,\n), E^ ^Skfik 



Nh 



(18) 



N/2 

HR = Y.''khVK 



Nh 



-N/2 



K = 



N/2- 



km 

sm I -— I a - 



N/2^l 



k=l 



N/2+k^N/2+j 



(26) 
(27) 
(28) 



Note two prominent facts: i) the single mode eigenenergies 
6:^ are the same for the L-chain and the i?-chain, because 
the two chains are identical, (ii) The Jordan Wigner oper- 
ators a'j of the i?-chain differ from the Jordan Wigner op- 
erators aj of the L-chain and total chain, because the lat- 
ter begin with spin 1, while the i?-chain begins with spin 
N/2-\- 1. We shall denote the eigenvectors and eigenvalues 
of as \r) = r2, . . . r7v/2) and Er, respectively: 



HR\r)=Er\r), 



N/2 



Nh 



(29) 



The states \n) are the Fock states associated to the fer- 
monic operators hk- For future reference we recall their 
properties 

h\hk\...nk ...) = n/e|...n/e ...) (19) 
Uk ...) = nfc(-l)^'=i' |... Uk - 1 ...) (20) 
hl\...nu...) = (l-nfc)(-l)^'- ".|...nfc + l...) (21) 



4 Initial quantum Entropy 

At times t < the entropy of the two non-interacting 
chains is given by the sum of their individual entropies. 
We proceed by calculating the quantum entropy of the 
left chain at inverse thermal energy /3. Since the two chains 
are identical this also gives the quantum entropy of right 
chain at the same inverse thermal energy p. 
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According to Eq. Q 
1 



m 



(30) 



Here L is the principal quantum number operator associ- 
ated to the L-chain. Its eigenvectors are the Fock states 
|/), and its eigenvalues are A^. The eigenvalues A/ are cal- 
culated in the following way. The energy eigenvalues Ei 
are ordered accordingly to their increasing values, so as to 
obtain a sequence 



Ei,<Ei,< - < El 



2N/2 



(31) 



where is the ground state, I/2) is the first excited 
state, . . . 1/2^/2) is the state of highest energy. Then A^^ = 
1, A^2 — 2, etc. 

The total initial quantum entropy 5*0, is given by: 



50 = 5^^+5^^. 



(32) 



5 Final quantum entropy 

Due to the assumption of sudden quench, at time t = 0+ 
the density matrix retains the initial form in Eq. ( 11 ). The 
final quantum entropy Si is therefore given by 



^1 



Xre-^^^^e-^^^^ In TV 



(33) 



Here TV is the principal number associated to the total 
chain, and i/n are its integer eigenvalues, obtained as de- 
scribed above for the smaller L-chain, with the difference 
that now one has to order the 2^ eigenvalues En. 

We next consider the basis \lr) = \l) ^ \r) formed by 
the direct product of the eigenbasis of and Hji: Using 
the resolution of the identity ^ \lr){lr\ = 1, the final 
quantum entropy reads: 

^ E.,^,.e-/^^^^e-/^^^^ In^^ P(n,/r) 



where 



Z{h)Z{pR) 



P{n,lr) = {n\lr){lr\n) . 



(35) 



In order to calculate P{n^lr) we consider the basis \s) = 
l^i, S2 . . . , sn) = |5i)(8)|52)(8)- • -^Isat) formcd by the direct 
product of the eigenstates \sj) the z component of each 



spin operator: (Jj\sj) = Sj\sj) with sj = ±. Employing the 
resolution of the identity twice we obtain: 

P{n,lr) = ^{n\s){s\lr){lr\s'){s'\n) 

= Y,{n\s){s'\n){lr\s'){s\lr). (36) 

The next crucial step in the calculation consists in ex- 
pressing the operators in terms of spin rising and 
lowering operators. For a single spin, say spin j, we have 
|+)(+| = a+a- , \-){-\ = a-a+ , |+)(-| = a+ , |+)(+| = 

(7 ■ which can be compactly written 



\s){s'\=5ss'<j'j<J-' + {l-5ss')cr' 



(37) 



where ^^5/ denotes the Kronecker symbol, and aj are 



given in Eq. ( |17[ ). Therefore: 

= Ws,s'y/cTp + (1 - 5,^,,X'^ . (38) 

j 

In order to calculate {n\s){s^\n) we express in terms 

of the fermionic operators bk of the total system. This can 
be accomplished by using the inverse sine transforms: 



an = 



N 



N- 



-^sin 

k=l 



kjiT 
TVTl 



and the following identities: 



(39) 

(40) 
(41) 



i-i 

a-^ = ]l{2a+a^ -l)a-\ 

where we set the convenient notations = apaj= aj. 
By plugging Eqs. (|39) [lo) [4l]) into Eq. (|38|, the operator 
is expressed in terms of the operators bk and 6^.. 
The wanted term (n\s ){s^\n) is calculated then by using 
the properties (Topl ) and the orthonormality condition 



In order to calculate {lr\s'){s\lr) we proceed by ex- 
pressing \s) as the direct product of two subchain states 
\s) = \si... sn) = |si . . . SN/2) (8) |sAr/2+i •••Sn) = \sl) 
\sii). Accordingly, the wanted term reduces to the product 
of two terms each pertaining to each subchain: {lr\s'){s\lr) 
= {l\s'L){sL\l)-{r\s'R){sR\r), The calculation of (/ 1 ) (sl | /) 
proceeds as the above calculation of {n\s){s'\n)^ with the 
only difference that now all quantities pertain to the smaller 
L-chain. Since the L and R chains are identical this auto- 
matically gives also the i?-chain term {r\s'j^){sR\r). 

In the Appendix we provide further details on how 
these calculations were implemented. 



6 Entropy change 

Using the method detailed above, we have calculated the 
change of quantum entropy 



AS = Si- So 



(42) 
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Fig. 2. (a): Entropy change, Eq. (42), as a function of Tl,Tr. Fig. 3. (a): Entropy change, Eq. (42), as a function of Tl,Tr 



(b): Entropy change, Eq. (42), as a function of Tr for various 
fixed Tl's. Here /i = 1, J = O.OOl/i and N = 8. 



(b): Entropy change, Eq. (42), as a function of Tr for various 
fixed Tl's. Here h = 1, J = h and N = 8. 



caused by the sudden quench, for different values of the pa- 
rameters defining our problem. These are the initial tem- 
peratures of the left and right chains, T^ n = Pj}^^ the 
length N of the total chain, and the interaction strength J. 
The field strength h was used as the unit of energy, and we 
adopted the convention that temperatures are measured 
in units of energy. In this unity we have /c^, Boltzmann's 
constant, equal to 1. 

Figure [2j top panel, shows a surface plot of the entropy 
change in a chain of length = 8 at J = O.Ol/i, as a 
function of Tl and Tr. Figure [2j bottom panel, shows for 
the same values of N = 8, J = O.Ol/i, the behavior of AS 
as function of Tr for various left chain temperatures Tl = 
0.5/i to 3.5/i. The quantum entropy change here is always 
positive and approaches a saturation value as Tr becomes 
very large. We see qualitatively different behaviors of AS 
as a function of Tr, depending on Tl- As Tl increases, AS 
changes from a monotonically increasing function of Tr to 
a monotonically decreasing function of Tli. The transition 
occurs for Tl of the order of the width of the spectrum of 
Hl (which takes on the value 2.6/i in this case). 

In Fig. |3] we report the same quantities as in Fig. [2] 
with the only difference that now J = h. For this value 
J = h the quantum entropy change, seen as a function of 



Tr features minima for small values of Tl and a mono- 
tonically decreasing behavior for large Tl. As with Fig. 
?? the transition occurs for Tl of the order of the width 
of the spectrum of Hl (which is 2h in this case). 

In both Fig. [2] and ^ we see a common feature. The 
values of AS corresponding to lower Tl is smaller than 
the value of AS corresponding to larger Tl at lower Tr, 
whereas it is vice- versa at larger values of Tr. We ob- 
served a similar behavior also at other values of J (e.g., 
for J = 0.0001,0.1,2,5) and AT = 2,4,6 for same value 
of all other parameters). In no case have we observed a 
negative change in quantum entropy. 



7 Thermalization 

After the quench, the system is in an out-of-equilibrium 
state. In order to quantify how far is the system from an 
equilibrium Gibbs state, one can employ one of the many 
metrics in the space of density matrices discussed in the 
literature ^28^. Among them the Hilbert- Schmidt distance: 

DHs[p,a] = VTr(p-(T)2 (43) 
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appears well suited to the problem at hand. The reason 
is that, in our problem, the Hilbert-Schmidt distance be- 
tween the after-quench density matrix pt and a Gibbs state 



(a) 



(44) 



where Zi{j3) = Tre~^^\ does not depend on time t. Fur- 
thermore it can be calculated by knowing the the initial 
density matrix po, Eq. (11 ), and the transition amplitudes 
P(n,/r), Eq. (36). That is, it can be obtained from the 
only knowledge of the (time independent) diagonal ele- 
ments 



{n\pt\n) = Y,P{n,lr) 



Ir 



^ — (SlEi ^ — (SfiEr 



(45) 



with no need to calculate the off-diagonal elements 
{n\pt\m)^ n ^ m. In fact: 



and 



n 



3-2/3£;„ 

WW' 



{ri\pt\n) , 



Trpl = Trpl = ^ 



^-2(3lEi ^-2pREr, 



(46) 

(47) 
(48) 

(49) 



Using Eqs. ( |43p9 ) we calculated the minimal distance 

Df^ = mmDHs[PuPp], (50) 

between the final state pt and the set of thermal Gibbs 
states. This gives both an estimate of how far the system 
is from equilibrium, and what the temperature T = 
is of the closest equilibrium, where ^ is the value of f3 for 
which the minimum distance D^^^^ is attained. 

Fig. [4] shows plots of T and D^^s (panels (a) and (b), 
respectively), as a function of Tl^Tr, at a low value of 
interaction strength J = 0.01 and chain length N = 
8. Panel (c) presents D^^^ as function of Tr for various 
fixed Tl's. Panel (b) indicates that better thermalization 
is achieved when Tl and Tr are closer. As can be seen 
from panel (a), when Tl = Tr it is T Tl = Tr. For 
fixed Tl (panel (c)) we observe that as Tr grows from 
zero, the minimal distance D^^^ first decreases, reaches 
then a minimum, and finally grows. The minimum is in 
correspondence to Tr ^ Tl as expected, and it is sharper 
for low Ti^'s. 

Fig. [5] is like Fig. [4] with the only difference that it is 
for a larger interaction strength J = h. As compared with 
Fig. |4] we observe here a different structure of the plot of 
^Ws ^ function of Tl and Tr, see panel (b). Within 
the range of Tl^ Tr considered in the plot, it appears that 
a smaller distance is achieved when either Tl^ Tr or both 
grow. This different structure is reflected in the curves 
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Fig. 4. (a): Temperature T of closest equilibri um , as a func 



tion o{Tl,Tr. (b): Minimal distance D'^f, Eq. (50 
tion of Tr.Tr. (c): Minimal distance D^^'^, Eq. ^ 
tion of Tr for various fixed Tl's. Here h — 1^ J — 
N = 8. 



) , as a func- 
I , as a func- 
O.OOl/i and 



presented in panel (c), presenting D^^^ as a function of 
Tr, for fixed T/^'s. The structure of the plot of T as a 
function of Tl^Tr^ panel (a), is qualitatively similar to 
the corresponding plot in Fig. [4j Note that for Tl = Tr, 
the temperature of the closest equilibrium is larger than 
Tl, of some amount of the order of h. This is because the 
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a function of Tr for various fixed Tl's. Here /i = 1, J = /i and 
AT = 8. 



quench injects an energy amount of the order J = h in 
the system. Similarly, in Fig. [Zj T was close to = Tr 
because the quench injects in that case an energy of the 
order J = O.Olh. 



We have investigated numerically the change in the quan- 
tum Hertz entropy of Eq. ([5|, caused by a quench of two 
spin-chains of different temperatures into a larger single 
chain. Although we cannot conclude that such changes are 
always positive, our numerics clearly suggests that this is 
the typical behavior, thus providing further support to the 
statement that the quantum Hertz entropy of Eq. ([5| is a 
proper quantum entropy for thermally isolated systems. 

We further quantified how far the system is from equi- 
librium after the quench, and estimated the corresponding 
temperature of the closest equilibrium. For those quenches 
ending very close to an equilibrium state, it becomes mean- 
ingful to assign the system the equilibrium temperature T, 
and the thermodynamic entropy Si. 
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Appendix 

As detailed in Sec.jsj the transition probabilities P(n, /r), 
involve the calculation of the expectation of the operators 
over the Fock states \n). Accordingly, we have de- 
tailed how these operators may be expressed in terms of 
the fermonic opera tors 6^, wh ose action on the Fock states 
is defined in Eqs. ( 19|2Q|21 ). In order to calculate those 
expectations we expressed the fermionic operators in 
matrix form. First we represented the Fock states \n) as 
tensorial product of single-spin states: 



|0) 



(51) 



For example, the Fock state |00) and |11) of a chain of 
N = 2 spins read 



|00) 



(52) 



and similarly for larger chains. In this representation, the 
searched fermionic operators are represented by the fol- 
lowing matrix tensorial products: 




1 terms 



k terms 



(55) 
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where cr± = (a^±a^)/2 are rising and lowering operators, 
expressed in terms of the Pauh matrices a^'^'^, and 1 is 
the 2x2 identity matrix. 

The calculation of P(n, Ir) further requires the calcu- 
lation of the expectations {1\sl){s']^\1) ^ {r\sR){s'^\r) . The 
calculation of these proceeds exactly in the same way de- 
tailed above, with the only difference that the chain length 
is now N/2 instead of N . With all these expectations one 
can calculate the probabilities P{n^lr)^ and, in turn, via 
Eq. p4| ), the final quantum entropy. 

The performance of the calculation can be greatly im- 
proved if one notices the following selection rules 



N N 

k=i j=i 

N/2 N/2 

k=i j=i 

N/2 N 

2^rfc7^ ^ Sj^N/2 =^ {r\sR)=0. 

k=l j=N/2-\-l 



(56) 
(57) 

(58) 



23. R. Dorner, J. Goold, C. Cormick, M. Paternostro, V. Ve- 
dral, Phys. Rev. Lett. 109, 160601 (2012) 

24. J. Gibbs, Elementary Principles in Statistical Mechanics 
(Yale U. P., New Haven, 1902), reprinted by Dover, New 
York, 1960. 

25. L. Landau, E. Lifschitz, Statistical Physics, 2nd edn. (Perg- 
amon, Oxford, 1969) 

26. E. Lieb, T. Schultz, D. Mattis, Ann. Phys. 16, 407 (1961) 

27. H.J. Mikeska, W. Pesch, Z. Phys. B 26, 351 (1977) 

28. J. Dajka, J. Luczka, P. Hanggi, Quantum Inf. Process. 10, 
85 (2011) 



Together with Eq. ( 36 ) these rules imply that the quench 
at time t = conserves the number of excitations. 
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